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Abstract 

We calculate twist-4 coefficient functions for the deep inelastic structure function F%(x, Q 2 ) 
associated to 4-gluon operator matrix elements for general values of the Bjorken variable 
x and study the numerical effect on the slope dF2(x,Q 2 )/dlogQ 2 . It is shown that these 
contributions diminish the strongly rising twist-2 terms towards small values of x. 
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1 Introduction 



The strong rise of deeply inelastic structure functions F 2 (x,Q 2 ) and Fl(x,Q 2 ) for small values 
of the Bjorken scaling variable x = Q 2 /2p.q is a consequence of both the structure of the twist-2 
(r = 2) evolution kernels at leading order and the behavior of the non-perturbative flavor singlet 
input distributions at scales Qq ~ 1 GeV 2 . Potentially this behavior violates unitarity which has 
to be guaranteed by higher order QCD corrections. As well-known, already the next-to-leading 
order corrections diminish the rising rate, see e.g. 00, however, they do not lead to a saturation 
°f F2,l(%, Q 2 ) as x — ■> 0. 

The large growth of F 2 (x, Q 2 ) in the double-logarithmic approximation with idealized 
5(1 — x)-like non-perturbative input distributions led Gribov, Levin and Ryskin [fj to propose a 
non-linear integro-differential equation for the gluon density, which obeyed saturation as x — > 0. 
In this scenario it is assumed that the non-perturbative n-particle gluon distribution is the 
n-th power of the single-particle gluon density. The evolution kernels are derived using the 
double-logarithmic approximation and an idealized factorization between the non-perturbative 
input densities and the evolution kernels is assumed. The negative sign of the non-linear gluon- 
recombination term is introduced referring to the Abramowskii-Gribov-Kancheli cutting rules Q 
from Regge theory. The resummation equation can be expanded into an infinite series of so- 
called fan-diagrams of 0(1/ (Q 2 ) k ). This equation is of the Glauber-type f| and not derived as a 
resummation of certain dominant classes of contributions emerging in the light cone expansion, 
although the term twist is used by some authors to denote the different l/(Q 2 ) k contributions.^ 

Later Mueller and Qiu || studied this recombination process too and derived a value for 
the recombination strength being yet unspecified in Ref. f|. Durand |7fl tried to reverse the 
mechanism leading to the leading order evolution equations, which, however, turned out to be 
not possible at this level.0 A first numerical illustration of the contributions due the resummation 
in || on the gluon density was pursued in Ref. §|. More detailed investigations, partly leaving 
the double-log level, were carried out in Refs. [10-12]. The non-linear resummation equation 
was solved in three different ways in [ITJ by showing existence and uniqueness of the solution of 
this equation. An exact value of the resummed gluon density in the saturation limit was derived 
analytically as well as the solution in the quasi-classic limit was studied leading to a derivation 
of a critical saturation line function of x and Q . 

In a forthcoming investigation of the resummation H by Bartels ||13|| , see also it was 
shown that even using the double-logarithmic approximation other than the fan-diagram con- 
tributions were present in the same order though suppressed by a power the number of colors, 
N c . Since the resummation is non-linear these perturbing terms being present in all orders spoil 
the resummation already in this approximation and quantitative predictions on their effect are 
hard to make. 

All these investigations have to be viewed in the general context of perturbative QCD. Here 
the question arises to which end the assumptions made can be validated in a rigorous approach. 
As well-known, the double-logarithmic approximation often yields much larger results than the 
complete calculation, cf. In a non-linear iteration quantitative effects due to this difference 
can be very large. To separate the non-perturbative input distributions and the perturbative 
evolution kernels and coefficient functions factorization has to be shown for the respective orders 
in l/(Q 2 ) k . In general a non-universal structure emerges both at the color and Lorentz-level 
comparing different orders in 1/Q 2 and the possibility to resum in general appears to be unlikely. 

2 Recent numerical estimates for NNL order Q continue this trend. 
3 More appropriately one should call these terms power-corrections. 

4 A similar mechanism is, however, applicable in the case of 2-nucleon reactions in a nucleus, 0. 
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Therefore, in a first step, the simpler question is asked for the structure of the 4 gluon — >• 
2 quark coefficient functions, which relates the 4-gluon operator matrix elements to the deep 
inelastic structure functions, cf. also Ref. jnj. Here a central question is that for the sign of 
these coefficient functions as a function of x, which in a sense was conjectured to be negative 
in the foregoing literature in the small x limit. It is clear that one neither can rely on the 
double-logarithmic approximation nor a small-x approximation in the sense of 'dominant poles', 
as well-known, see Refs. [[17], [L^, [18] , due to the interplay of small x and large x contributions via 
the Mellin-type integrals between the coefficient functions and the respective operator matrix 
elements. Instead one has to perform a calculation for the whole range of x. 

In this note we present numerical results on the effect of main contributions to the twist-4 co- 
efficient function driven by a two-particle gluon distribution G 2 for the slope dF 2 (x, Q 2 )/d\ogQ 2 . 
This quantity is of interest since one may generally expect that higher twist effects are more 
easily detectable in this distribution than in the integrated structure function F 2 (x, Q 2 ). After 
discussing the structure of the non-perturbative distribution G 2 in section 2 the coefficient func- 
tion is derived (section 3). Numerical results are presented in section 4 and section 5 contains 
the conclusions. 



2 4— Gluon Operator Matrix Elements 

The non-perturbative input distributions for a 4-gluon state being linked to a quark box contains 
four color, four vector indices, and three scaling variables in general. A priori nothing is known 
about this quantity except the weight factor Q\jQ 2 relative to the twist 2 contributions defined 
by the non-perturbative parameter Q\ = 1/R 2 . Here R denotes a twist-4 screening length. 
In analogy to the single particle gluon density one might define the x-dependent part for a 
multi-gluon density depending on n momentum fractions 

G {n) (x 1 ,...,x n )^A n x^...x^{l-x 1 -...-x n ) (1) 



similar to Ref. pi. One may further assume a = aij,Vi for simplicity. If one attempts to relate 



the four-gluon distribution 

G 2 (xi, x 2 ; x[, x' 2 ) = G 2 (xi, x 2 \ x' 1 ,x' 2 )5(xi + x 2 - x[ - x 2 ) (2) 

to a product of two single-gluon distributions Gi(xj) it becomes clear that this is only possible 
by discarding one of the x/s or relating it to the other two. In a first attempt one still might 
want to work in such an approximation and even assume that G 2 is found as the simple product 
of two single gluon densities setting the unknown non-perturbative correlator x{ x ii x 2i x z) t° 1 
to obtain first numerical results. 

It is evident that on the level of twist-4 the general description requires the introduction of 
several 4-particle gluon distributions, according to the different color and Lorentz structures, 
which are not related. This further complicates the determination of these quantities from ex- 
perimental data. The permutation behavior of the arguments of G 2 being linked to a certain 
piece of the coefficient function depends on the transformation properties of the latter concern- 
ing both color and the momentum fractions X{. Given the fact that the correlators \ are yet 
unknown, it is not even clear a priori, whether G 2 (xi,x 2 ,x^) is positive everywhere. With all 
these reservations in mind we are going to follow a simplified approach in the present paper. We 
will assume X\ = x[ = x 2 = x' 2 and x = 1- This leads to 

GtM = ^G?(xO (3) 



2 



which guarantees the positivity of G<i- Since the ansatz (|^) is not of general validity one has 
to restrict the investigation using special color and Lorentz contractions in the factorization of 
the coefficient function and the hadronic matrix elements below. Following || we use for the 
representation of the two-particle (momentum) gluon density 



G a (rci) = %Q\x\G\{x x ) . 



(4) 



3 Coefficient Function 



The twist-4 coefficient function describing the 4-gluon — > 2-quark amplitude is calculated to 
dp 2 j_/p\ accuracy using time-ordered perturbation theory [OT. A first investigation of the prob- 
lem has been performed in [^l], p2| . The contributing terms are illustrated in Figure 1 (direct 
terms) and Figure 2 (interference terms). Here the grey ovals symbolize the set of diagrams 
being shown in Figure 3. The transverse momentum p± occurs in the internal loop. 

t -> 




Figure 1: Direct diagrams contributing to (|5|). The grey oval symbolizes the set of diagrams in Figure 3. 
Orthogonal dashed lines stand for the time ordering. The separated white ovals symbolize the two parts 
of the non-perturbative 4-gluon density. x\,X2,x[ and x' 2 are the longitudinal momentum fractions, with 



X\ + x 2 



0. The circle stands for the forward subprocess 7* + q — > 7* + q through which the 



virtual photon couples to the amplitude. 
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Figure 2: Interference diagrams associated to the process in Figure 1. 

The white ovals in Figures 1,2 stand for the respective parts of the non-perturbative distribution 
G2, which are collinearly factorized. In the calculation we assumed x\ = x[ and X2 = x' 2 , which 
is the natural choice, if one aims on using the ansatz Ga^Ei,^) = xiGi(x\)x2Gi(x2)- Finally 
X\ and X2 were, furthermore, set equal, in a first approximation, according to the choice Eq. (Q) 
for G 2 . 
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Figure 3: Diagrams symbolized by grey ovals in Figure 1 and 2. 



We have independently recalculated the coefficient functions for these processes. For the 2-^2 
coefficient function one obtains : 



a 



4^2;2->2/ 
G^q 



[X\ } X\, X^, X]] X 



a 



4^2;2->2, 
G^q '■ 



1 (2xi 



X 



96 



x\ 



[Ux 2 - 3xx 1 + 18x*j 



(5) 



Under the above assumptions the 3rd diagram in Figure 3 does not contribute. Eq. (|5|) can be 
given a simple form introducing the variable z = xjx\. Note that the Mellin poles in this variable 
are situated at iV = 0, —1, . . . , —4 and are of similar strength. Non of these poles is therefore 
dominant. For the direct process the kinematic range of the variable X\ is x\ e [x/2, 1/2]. 

The interference contributions, Figure 2, consist out of the coefficient function for the 1 — > 3 
and 3 — > 1 process and the factorized non-perturbative distributions. One may show, cf. ET 



that these coefficient functions are related to that of the 2^2 process by a sign-change in the 
whole x range. 



a 



4^2;2^2 
G^q 



-a 



4^2;1^3 
G^q 



{x\i X\] X-y, X±] x) 



-a 



4^2;3^1 
G^q 



(Xi 

)• (6) 



However, the kinematic range of X\ is now x\ e [x, 1/2]. Arguments were given in Ref. [23| that 
one may choose the same non-perturbative input for the direct and interference terms. This 
hypothesis is being adopted here. For the slope dF 2 (x, Q 2 )/d\ogQ 2 one finally obtains : 



dF 2 (x,Q 2 ) 
d\og(Q 2 ) 



dFl= 2 \ 



x,Q 2 



d\og(Q 2 ) 



+ 



2nJ Q 2 j x /2 



1/2 



X 



f'-' j ( C^q 2 ~* 2 (x 1 ,x)G 2 (x 1 ) 



1/2 



'G^q 



x 1 ,x)G 2 (x 1 ) , (7) 



where 



dF 2 T=2 \ 



x,Q 2 



<91og(Q 2 ) 

and eg) denotes the Mellin convolution 



[x 



4 Numerical Results 

Numerical values for the slope dF 2 (x, Q 2 ) / d\ogQ 2 are given in Figure 4. Here we compare the 
results in leading order QCD for the twist-2 contribution to those obtained including the twist-4 
term Eq. (|7]). For the twist-2 contribution and the non-perturbative distribution G2{z,Q 2 ) we 
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refer to the parameterization [24] and Eq. (|), respectively. Figure 4a and b depict the results for 
Q 2 = 5 GeV 2 and 10 GeV 2 , choosing different screening lengths R. Under the above assumptions 
the twist -4 correction is negative in the small-x range, i.e. the screening term in Eq. fl7|) wins 
against the anti-screening term. The latter one leads to a very small positive contribution only 
at large values of x. The effect shrinks with rising values of Q 2 , and becomes visible only for x 
values below x ~ 10~ 4 at perturbative scales of Q 2 ~ 4 GeV 2 . Therefore this effect cannot be 
seen yet in the kinematic range of HERA. Although one might find arguments for the size of the 
screening parameter R, it is a non-perturbative quantity for which a general prediction within 
QCD is still missing. We gave illustrations for R 2 = 5 and 2 GeV 2 . In the latter case a clear 
effect is visible for x < 2 • 10" 6 and Q 2 ~ 5 GeV 2 . 

The coefficient function Eqs. (|5|,§) applies to the whole 2-range with z = xjx\. One might 
try to find an effective small-x approximation. The function is a polynomial of degree 4 in 
the variable z. Approximating it successively up to its z°,z,z 2 . . . -contribution one may study 
the effect of the individual terms. Figure 5 shows, that a sufficient representation is only given 
taking all contributions into account since all approximants diverge for x — > 0, i. e. an effective 
small- x approximation to this function does not exist. 

5 Conclusions 

A study has been performed to the slope of F 2 due to twist-4 coefficient functions for the 
process AG 2q. The present calculation was performed in time-ordered perturbation theory 
and limited to dp\/p\ accuracy of the respective Feynman diagrams. Numerical results were 
choosing a special ansatz for the non-perturbative two-particle gluon density G2. The new 
contributions contain anti-screening and screening terms, the latter of which dominate for smaller 
values of x and diminish the growth of the slope due to the twist-2 contributions. The numerical 
results show that this effect is of significant size only below x ~ 2 • 10~ 6 for Q 2 ~ 5 GeV 2 . This 
region is yet beyond the kinematic range which can be probed at the ep collider HERA but may 
be investigated at future lepton-hadron facilities operating at larger cms energies. 
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Figure 4b: Same as figure la for Q 2 = lOGeV 2 . 



8 




Figure 5: Comparison of the slope dF2(x,Q 2 )/d\ogQ 2 at Q 2 = 5GeV 2 and twist-4 mass scale 
R 2 = 2 GeV~ 2 , Eq. (0) (full line) with the corresponding results obtained approximating the coefficient 
function Eq. ([5]) by the sequence of contributing powers. 1st: z°, 2nd: z etc. (dotted lines). Dash- 
dotted line: twist-2 contribution. 



9 



